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[57] ABSTRACT 

The present invention is a system for chemometric analysis 
for the extraction of the individual component fluorescence 
spectra and fluorescence lifetimes from a target mixture. The 
present invention combines a processor with an apparatus 
for generating an excitation signal to transmit at a target 
mixture and an apparatus for detecting the emitted signal 
from the target mixture. The present invention extracts the 
individual fluorescence spectrum and fluorescence lifetime 
measurements from the frequency and wavelength data 
acquired from the emitted signal. The present invention uses 
an iterative solution that first requires the initialization of 
several decision variables and the initial approximation 
determinations of intermediate matrices. The iterative solu- 
tion compares the decision variables for convergence to see 
if further approximation determinations are necessary. If the 
solution converges, the present invention then determines 
the reduced best fit error for the analysis of the individual 
fluorescence lifetime and the fluorescence spectrum before 
extracting the individual fluorescence lifetime and fluores- 
cence spectrum from the emitted signal of the target mixture. 

16 Claims, 6 Drawing Sheets 
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CHEMOMETRIC ANALYSIS FOR 
EXTRACTION OF INDIVIDUAL 
FLUORESCENCE SPECTRUM AND 
LIFETIMES FROM A TARGET MIXTURE 

5 

The U.S. Government has a paid-up license to certain 
technologies disclosed in this invention and the right in 
limited circumstances to require the patent owner to license 
others on reasonable terms as provided for by the terms of 
the following contracts: NAS1-20426 and NAS1-20162 10 
awarded by NASA; DAAH01-91-R198 awarded by US 
Army Missile Command; DAAL06-92-C-0014 and 
DAAD07-91-C-0127 awarded by US Army White Sands 
Missile Range; DAAA15-93-C-0034 and DAAM01-94-C- 15 
0033 awarded by US Army Chemical and Biological 
Defense Agency; and F41624-95-C-6010 and F41624-97- 
C-6029 awarded by US Air Force Armstrong Laboratory. 

BACKGROUND OF THE INVENTION 20 

1. Field of the Invention 

The present invention generally relates to fluorescence 
measurements. More specifically, the present invention 
relates to frequency domain measurements of the fluores- 2 s 
cence lifetime and the fluorescence spectrum. 

2. Description of the Related Art 

The process of fluorescence occurs when a substance, 
such as a molecule, absorbs light at one wavelength (or 
energy), and then emits light at a longer wavelength (or 
lower energy). A slight time delay occurs from when the 
substance absorbs light and when the substance re -emits 
light at the longer wavelength. This time delay is known as 
the fluorescence lifetime. 35 

FIG. 1 depicts the fluorescence process schematically for 
a diatomic molecule, which is a molecule that is comprised 
of two atoms. There are several “types” of energy that are 
internal to a diatomic molecule, one of which is vibrational 

40 

energy. The classic example of vibrational energy is imag- 
ining the two atoms connected together by a spring where 
the atoms oscillate back and forth along the axis of the 
spring. As the two atoms approach each other, they experi- 
ence repulsion due to the proximity of their negatively 45 
charged electron clouds. As the atoms pull apart from one 
another, they experience an attraction that one can imagine 
as the result of the attractive forces between the positively 
charged nucleus of one atom for the negatively charged 
electron cloud of the other atom, and vice versa. The 50 
distance between the two atoms is known as the internuclear 
distance, which constantly changes as the two atoms oscil- 
late back and forth. As the internuclear distance decreases, 
a rise in potential energy occurs due to the repulsive forces; 55 
as the internuclear distance increases, a rise in potential 
energy again occurs due to the attractive forces. FIG. 1 
illustrates this difference in potential energy by line 11 in the 
potential energy “well” 10 and by line 13 in the potential 
energy “well” 12 where line 11 and 13 plot the potential 60 
energy as a function of internuclear distance. 

A group of horizontal lines 20, 22, and 24 appear in each 
of the potential energy wells 10 and 12. Each horizontal line 
represents the individual vibrational energy states possible 65 
for the diatomic molecule. Quantum mechanics requires that 
the vibrational frequency of the spring be within certain 


2 

values. In other words, the spring may oscillate at frequency 
V, and the spring may oscillate at frequency ‘b’, but it is a 
physical impossibility for the spring to oscillate at any 
frequency between ‘a’ and ‘ b \ In FIG. 1, any two adjacent 
horizontal lines would represent ‘a’ and ‘ b \ 

The different potential energy wells 10 and 12 represent 
another type of energy in the diatomic molecule: electronic 
energy. Different electronic energies occur when the mol- 
ecule absorbs energy in such a fashion that it causes an 
electron to move to a higher energy configuration within the 
molecule. The classical example of the energy absorption is 
the 'changing the spring’ that connects the two atoms. If we 
“add” electronic energy to the potential energy well 10, the 
energy raises the potential energy well 10 to the energy level 
of potential energy well 12. 

With this background, we can use FIG. 1 to describe the 
process of fluorescence and fluorescence lifetime for a 
diatomic molecule. The example molecule is originally in a 
state depicted by the lower potential energy well 10. The 
molecule absorbs a photon of the correct amount of energy 
that induces an electron to move to a higher electronic 
energy state, represented by the upper potential energy well 
12. The diatomic molecule also tends to undergo a change 
from a Tower vibrational state’ (the horizontal lines 20 and 
22) in the lower potential well 10 to a higher vibrational state 
in the upper potential well 12. This initial excitation of the 
molecule is shown by line 14. Through any one of numerous 
possible processes, the higher vibrational state 24 in the 
upper well decays to a lower vibrational state 22 in the upper 
well (as shown by line 18). The process of non-radiative 
decay to a lower vibrational energy level in the upper well 
occurs very rapidly. After the non-radiative vibrational 
relaxation, the electron will want to return or revert back to 
a lower energy state as shown by line 16 (in other words, the 
electron wants to move back to where it was before the 
whole process started), and in so doing, the molecule emits 
a photon of light (with a lesser amount of energy than the 
excitation photon). This lowering of the energy state and the 
emission of a photon is the process of fluorescence. The 
process of emitting the photon has an associated time delay, 
which is the fluorescence lifetime. 

In practice, we typically probe many different types of 
molecules at once with the excitation light pulse. FIG. 2 
illustrates the case where a single short excitation pulse of 
light 30 is absorbed by a sample of identical molecules all 
at once. The fluorescence decay curve 32 resulting from a 
typical fluorescence response of a sample of identical mol- 
ecules is exponential in nature because not all of the iden- 
tical molecules emits its fluorescence photon at precisely the 
same time. The exponential decay follows a mathematical 
function so that we can calculate the fluorescence lifetime, 
t: 

where t 0 is the time the excitation pulse, I 0 is the initial 
fluorescence and I(t) is the observed fluorescence intensity 
as a function of time. 

In the cases with multiple molecular types (where more 
than one lifetime decay is present within a sample of a target 
mixture), the exponential decay seen in FIG. 2 would appear 
as a sum of exponential functions. Prior art fluorometers 
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typically use a type of time -correlated (or time -resolved) 
system that count the emitted photons (from an excited 
molecule) in order to measure the fluorescence lifetime. 
Other prior art systems add the ability to take the fluores- 
cence lifetime measurements in the frequency domain by 
modulating the dynode of a photomultiplier tube, followed 
by a mixing and correlation procedure. These prior art 
systems are cumbersome, time consuming, and complicated 
to operate. The present invention overcomes the limitations 
of the prior art systems by utilizing a novel technique to 
measure the fluorescence lifetime and spectrum. Instead of 
irradiating the target sample with a single short pulse of light 
(photon counting), the present invention continuously irra- 
diates the target sample with a light source whose amplitude 
modulation frequency is stepped with time. This technique 
allows us to use a chemometric analysis to automatically 
extract the lifetimes from the 'phase delay’ and 'intensity vs. 
frequency’ characteristics of the emitted light. 

SUMMARY OF THE INVENTION 

The present invention is a system for chemometric analy- 
sis for the extraction of the individual fluorescence spectrum 
and fluorescence lifetime from a target mixture. The present 
invention comprises a processor with an apparatus for 
generating an excitation signal to transmit at a target mixture 
and an apparatus for detecting the emitted signal from the 
target mixture. One embodiment of the present invention 
uses a processor that comprises a computer that extracts the 
individual fluorescence spectrum and lifetime measurements 
from the frequency and wavelength data acquired from the 
emitted signal. The present invention first determines the G 
and S matrices from the frequency and wavelength vectors. 
A renormalization of w -1 S x to G x occurs next due to the 
differences in units. Next, the present invention determines 
the initial U, V, g and w -1 s approximations. The present 
invention uses an iterative solution that first requires the 
initialization of the decision variables including r\n. The 
iterative solution compares the decision variables for con- 
vergence to see if further approximation determinations of 
the U, V, g and w -1 s are necessary. When the solution 
converges, the present invention then determines the 
reduced best fit error for the analysis of the individual 
fluorescence lifetime and the fluorescence spectrum. And 
finally, the system of the present invention extracts indi- 
vidual component fluorescence lifetimes and fluorescence 
spectra from the emitted signal of the target mixture. 

The present invention additionally includes a method and 
apparatus for generating and mixing signals for frequency- 
domain lifetime and spectral fluorometry. The present inven- 
tion comprises a signal generator that generates a driving/ 
reference signal that modulates the amplitude and/or the 
frequency of the driving/reference signal over time and a 
signal generator that generates the mixing signal that modu- 
lates the amplitude and/or the frequency of the mixing signal 
over time. The driving/reference signal generator drives an 
excitation signal generator that generates the excitation 
signal. The present invention then directs and transmits the 
excitation signal at a target mixture, which absorbs the 
energy from the excitation signal. The property of fluores- 
cence causes the target mixture to emit an emitted signal that 
the present invention detects with a signal detector. The 
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present invention uses a mixer that mixes the emitted signal 
with the driving/reference signal to produce a processor 
reference signal. Another mixer mixes the emitted signal 
with the mixing signal to produce a data signal. The present 
5 invention then uses a processor to compare the processor 
reference signal with the data signal by analyzing the 
differences in the phase and the differences in the amplitude 
between the two signals. The processor then extracts the 
fluorescence lifetime and fluorescence spectrum of the emit- 
ted signal from the phase and amplitude information using 
a chemometric analysis. 

DESCRIPTION OF THE DRAWINGS 

15 To further aid in understanding the invention, the attached 
drawings help illustrate specific features of the invention and 
the following is a brief description of the attached drawings: 

FIG. 1 illustrates the fluorescence process for a diatomic 
20 molecule. 

FIG. 2 illustrates the absorption and emission of energy 
for the fluorescence process. 

FIG. 3 illustrates the chemometric extraction process of 
25 the present invention for fluorescence spectrums and life- 
times. 

FIG. 4 is a block diagram of the present invention. 

FIG. 5 is a block diagram of a component of the present 
30 invention for generating an excitation signal. 

FIG. 6 is a block diagram of an embodiment of the present 
invention for processing an emitted signal. 

FIG. 7 is a block diagram of the present invention that 
35 illustrates the electrical pathways of the excitation signal and 
the emitted signal. 

FIG. 8 is a block diagram of an embodiment of the present 
invention for a chemometric analysis of the emitted signal to 
determine the fluorescence spectrum and fluorescence life- 
40 times. 

FIG. 9 is another block diagram of an embodiment of the 
present invention for a chemometric analysis of the emitted 
signal to determine the fluorescence spectrum and fluores- 
45 cence lifetimes. 

FIG. 10 illustrates the present invention’s ability to dif- 
ferentiate spectra in a target mixture. 

DETAILED DESCRIPTION OF THE 
50 INVENTION 

The present invention is a system for chemometric analy- 
sis for the extraction of the individual component fluores- 
cence spectra and fluorescence lifetimes from a target mix- 
55 ture. Additionally, the present invention is a method and 
apparatus for generating and mixing signals for frequency- 
domain lifetime and spectral fluorometry. This disclosure 
describes numerous specific details that include specific 
processes, structures, and circuits in order to provide a 
60 thorough understanding of the present invention. For 
example, the present invention describes the extraction of 
the fluorescence spectra and fluorescence lifetimes from a 
target mixture. However, the practice of the present inven- 
65 tion includes the extraction or analysis of other physical 
properties of a target mixture other than the previously 
described ones. One skilled in the art will appreciate that one 
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may practice the present invention without these specific 
details. Additionally, this disclosure does not describe some 
well known processes such as Fourier transforms, 
eigenequations, eigenfunctions, eigenvalues, transforms, or 
best fit analyses in detail in order not to obscure the present 
invention. 

The present invention utilizes a novel system that uses a 
chemometric analysis to determine the fluorescence life- 
times and the fluorescence spectra of a target mixture. 
Instead of irradiating the target sample with a single short 
pulse of light (photon counting) as other prior art systems, 
the present invention continuously irradiates the target mix- 
ture with a light source whose amplitude modulation fre- 
quency is stepped with time. FIG. 3 illustrates the process of 
the present invention. The excitation light or signal 42, 
which is the signal the present invention uses to continu- 
ously irradiate the target mixture, is amplitude and phase 
modulated in a sinusoidal fashion. As previously discussed, 
the emitted signal 40 is the photon emitted (emission 
intensity) from the molecule as a result of the fluorescence 
process. The emitted signal 40 will appear as a delayed 
sinusoidal-intensity light of a longer wavelength. Embedded 
within the emitted signal is the information necessary to 
extract the fluorescence lifetime and spectrum. However, 
one must use a chemometric analysis to extract the fluores- 
cence lifetimes from the 'phase delay <D’ of the emitted light 
(compared to excitation signal) and the 'intensity vs. wave- 
length’ (or spectrum) characteristics of the emitted light 
(compared to the excitation signal). 

FIG. 4 is a block diagram of the apparatus component of 
the present invention. A processor 50 couples to a data 
acquiring component 52 through a data input/output path 68. 
The preferred embodiment of the processor 50 comprises a 
computer. The data acquiring component comprises a gen- 
erator component 54 and a receiver component 56. Coupled 
to the generator component 54 is an excitation signal gen- 
erator 58. The excitation signal generator produces an exci- 
tation signal (or light or photon) that is transmitted to a target 
mixture 66. The target mixture 66 absorbs the electronic 
energy of the excitation signal as the molecules of the target 
mixture move to a higher electronic energy state. During the 
fluorescence process, the molecules of the target mixture 
emit photons (an emitted signal 64) as the molecules move 
to a lower electronic energy state. A detector 60 detects the 
emitted signal 64 from the target mixture 66. Coupled to the 
detector is the receiver component 56, which processes the 
returned emitted signal 64 into a reference signal and a data 
signal for the processor 50. The processor 50 then uses a 
chemometric analysis to determine the individual fluores- 
cence spectrum and fluorescence lifetime of one or more 
fluorophores within the target mixture 66. 

FIG. 5 is a block diagram of the generator component 54 
of the present invention for generating the excitation signal. 
The present invention uses heterodyning techniques to pro- 
duce two sinusoidal RF signals, a driving/reference signal 
90 and a mixing signal 92. The present invention modulates 
the frequency of the two signals from 10 MHz to 200 MHz. 
One skilled in the art will appreciate that the present 
invention could vary the signals over a much larger fre- 
quency range. The preferred embodiment of the present 
invention generates the two signals with a frequency differ- 
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ence of 10 kHz. Another embodiment of the present inven- 
tion generates the two signals using an adjustable offset 
frequency where the offset frequency is set through the 
present invention’s control software. 

5 Referring back to FIG. 5, the present invention generates 
the driving/reference signal with a driving/reference signal 
generator 71 and generates the mixing signal with a mixing 
signal generator 73. The driving/reference signal generator 
71 comprises an oscillator 74, a phase-locked-loop (PLL) 
70, a mixer 78, and a signal splitter 72. The mixing signal 
generator 73 comprises an oscillator 76, the phase-locked- 
loop (PLL) 70, a mixer 80, and the signal splitter 72. The 
present invention generates the two signals by mixing a 
15 signal from the PLL 70 with the signals from the two fixed 
oscillators 74 and 76. The preferred embodiment of the 
present invention operates oscillator 74 at 290 MHz and 
oscillator 76 at 290.01 MHz. One skilled in the art will 
appreciate that other operating frequencies for oscillators are 
possible. The preferred embodiment of the present invention 
can set the output of the PLL 70 to a frequency ranging from 
300 to 560 MHz, which the present invention controls 
through the control software of the present invention 
25 through the digital I/O datapath 68. 

Coupled to the driving/reference signal generator 71 is a 
filter 82 and an amp 86. The filter 82 filters out any unwanted 
components from the driving/reference signal, and the amp 
86 provides additional gain to the signal if necessary. 
Coupled to the mixing signal generator 73 is a filter 84 and 
an amp 88. The filter 82 filters out any unwanted compo- 
nents from the mixing signal, and the amp 86 provides 
additional gain to the signal if necessary. The present 
35 invention uses the driving reference signal 90 as a phase and 
amplitude reference signal on the receiver component 56 (of 
FIG. 6). And, the present invention sends the mixing signal 
92 to the receiver component 56 where it is split. 

The present invention uses the driving/reference signal 90 
to directly drive and modulate an excitation signal generator 
58 (of FIG. 7) that may comprise a laser diode, a light 
emitting diode (LED), or a Deuterium lamp. The present 
invention can also use the modulation to directly modulate 
45 an electro -optic crystal for use in external modulation of gas 
laser sources or other continuous emitting excitation light 
sources. 

FIG. 6 is block diagram of the receiver component 56 of 
the present invention for processing the returned emitted 
signal. The present invention receives the modulated emitted 
signal 64 from the target mixture through a detector (60 of 
FIG. 7). The receiver component 56 receives the detected 
emitted signal 96 from the detector through an amplifier 97 
55 and receives the driving/reference signal 90 through an 
amplifier 99. The receiver component 56 additionally 
receives the mixing signal 92 through a splitter 98. The 
present invention then mixes the detected emitted signal 96 
with the mixing signal 92 with a mixer 102 and mixes the 
60 driving/reference signal 90 with the mixing signal 92 with a 
mixer 100. The resulting signals contain the sum and dif- 
ference frequencies of the mixer’s (100 and 102) input 
signals. The present invention generates a data signal 114 at 
65 10 kHz by filtering the mixed detected emitted signal 
through a filter 106 with an amplifier 110 to eliminate the 
high frequency components. The present invention addition- 
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ally generates a processor reference signal 112 at 10 kHz by 
filtering the mixed driving/reference signal through a filter 
104 with an amplifier 108 to eliminate the high frequency 
components. One skilled in the art after reading this disclo- 
sure will appreciate that one can use frequencies other than 
10 kHz to process signals. An analog to digital converter 118 
converts the data signal 114 into digital data for transmittal 
over the digital I/O path 68 to the processor (50 of FIG. 4 ), 
and an analog to digital converter 116 converts the processor 
reference signal 112 into digital data for transmittal over the 
digital I/O path 68 to the processor. The present invention 
then processes the two input signals to determine their 
relative amplitudes and phases because the amplitude and 
phase of the 10 kHz signals are directly proportional to the 
amplitude and phase of the high frequency modulation 
signals. As the modulation frequency is stepped from 10 
MHz to 200 MHz, the processor reference signal 112 and the 
data signal 114 remain at 10 kHz. 

FIG. 7 is a block diagram of the present invention that 
illustrates the electrical pathways of the excitation signal and 
the emitted signal. The present invention generates the 
driving/reference signal with the driving/reference signal 
generator 71 and generates the mixing signal with a mixing 
signal generator 73 . The present invention uses the driving/ 
reference signal 90 to directly drive and modulate an exci- 
tation signal generator 58 that may comprise a laser diode, 
a light emitting diode (LED), or a Deuterium lamp. The 
present invention can also use the modulation to directly 
modulate an electro-optic crystal for use in external modu- 
lation of gas laser sources or other continuous emitting 
excitation light sources. The excitation signal generator 58 
produces an excitation signal 62 that is transmitted to a 
target mixture 66 . A photodetector 120 detects the excitation 
signal 62 on its way to the target and regenerates the 
driving/reference signal 90 for retransmission to the receiv- 
ing component 56 of the present invention. The target 
mixture 66 absorbs the electronic energy of the excitation 
signal as the molecules of the target mixture move to a 
higher electronic energy state. During the fluorescence 
process, the molecules of the target mixture emit the emitted 
signal 64 as the molecules move to a lower electronic energy 
state. A detector 60 detects the emitted signal 64 from the 
target mixture 66 . The detector 60 could comprise a photo - 
detector. Additionally, wavelength dispersion of the emitted 
signal can be accomplished by using a spectrophotometer, 
monochromator or low-pass, high-pass and notch optical 
filters. From the detector 60 , the present invention mixes the 
signals as previously described for FIG. 6 above. The 
present invention then uses a chemometric analysis to deter- 
mine the individual fluorescence spectrum and fluorescence 
lifetime of one or more fluorophores within the target 
mixture 66. 

The present invention comprises several techniques or 
alternative embodiments to “target” the target mixture. One 
embodiment of the present invention uses a target compart- 
ment for holding the target mixture within a chamber. A 
container may contain the target mixture, or the mixture may 
be free floating, in a fluid state, within the chamber. Another 
embodiment of the present invention uses an optical probe 
that combines the signal generator 58 and the detector 60 
into a single probe. A user of the present invention would 
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then probe or place the optical probe into the target mixture. 
And, another embodiment of the present invention may use 
a laser to target a free floating target mixture, such as the 
emission from a car (at a distance). One skilled in the art will 
5 appreciate that other embodiments of targeting the target 
mixture are possible after reviewing the disclosure of the 
present invention. 

FIG. 8 is a block diagram of an embodiment of the present 
10 invention for a chemometric analysis of the emitted signal to 
determine the fluorescence spectrum and fluorescence life- 
times. Before discussing the block diagram of FIG. 8, we 
first need to discuss the mathematical foundation used for 
the chemometric analysis. Fluorescence occurs when a 
molecule, excited by a photon to a higher electronic state, 
decays non-radiatively to a lower vibrational level of the 
excited state, and then decays radiatively to a vibrational 
level of its ground electronic state. The non-radiative decay 
20 is assumed to be ‘fast’. The characterization of the radiative 
decay is the fluorescence lifetime, t. And as previously 
stated, the wavelength of the fluorescent photon is greater 
(meaning less energy) than that of the exciting photon. 

25 If we apply a ‘monochromatic’, short light pulse to the 
molecules of a fluorescent material at a time t 0 , ideally 
represented by: 

E(i)=E 0 -b(t-Q (l) 

30 where E 0 is the total energy in the pulse and E(t) is the 
intensity of the incident pulse as a function of time, we find 
that the fluorescent intensity at a fixed wavelength has the 
form: 

35 (2) 

where 

I 0 =aE 0 (3) 

40 and a has appropriate units. Thus, by measuring the decay 
of the fluorescent intensity with time and fitting a straight 
line (linear least squares fit) to a graph of log(I(t)) vs time, 
we can determine the lifetime of a fluorophore. However, if 
more than one fluorophore is present, the intensity variation 
45 with time depends on the concentrations, c„, and lifetimes, 
t„, of the individual fluorophores: 

I(t ) = / 0 ^ c n e- [t - f o )/T » (4) 

n 

50 

Now the intensity decay curve must be fit by a non-linear 
least squares process, which requires initial guesses for the 
number of fluorophores, r, their concentrations, and their 
55 lifetimes. In general, convergence of the solution cannot be 
guaranteed as it depends on the accuracy of the initial 
guesses. 

Weber showed that decay constants and relative intensi- 
ties of an arbitrary number of independent components in a 
60 heterogeneous fluorescence emission could, in principle, be 
determined by measuring phase shift and relative modula- 
tion of the total fluorescence at an equal number of frequen- 
cies. See Weber, G., Resolution of the Fluorescence Life- 
65 times in a Heterogeneous System by Phase and Modulation 
Measurements, J. Phys. Chem., 85, 949-953 (1981). At each 
frequency, the phase and modulation measurements yield 
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the complex Fourier transform of the fluorescence impulse 
response, G+iS. The moments of a distribution of lifetimes 
are found as linear combinations of the G’s and S’s and 
Prony’s method is used to obtain lifetimes and concentra- 
tions of the components from the moments. Due to its 5 
sensitivity to noise, this method is really not practicable for 
more than binary mixtures. 

Lackowicz and his coworkers developed an analysis of 
frequency domain phase and modulation data based on a 10 
weighted non-linear least square analysis. See Lackowicz, J. 

R., et al., Analysis of Fluorescence Decay Kinetics from 
Variable-Frequency Phase Shift and Modulation Data, Bio- 
phys. J., 46, 463-477 (1984). By taking data at several 
emission wavelengths, they were able to resolve binary 15 
mixtures whose components differed in lifetime by 30% and 
ternary mixtures the ratio of whose range of lifetimes 
spanned a decade. Best fit is determined by minimizing a 
reduced chi-squared function. Poor initial guesses can still 20 
produce unpredictable results. 

Working with an Excitation-Emission-Frequency Array, 
McGown and her co workers employed a principal compo- 
nents analysis followed by a non-linear least squares mini- 
mization on an over determined set of data to extract the 25 
spectra and lifetimes of components in a mixture. See 
Burdick, D. S., et al., Resolution of Multicomponent Fluo- 
rescent Mixtures by Analysis of the Excitation-Emission- 
Frequency Array, J. Chemometrics, 4, 15-28,(1990). In 3Q 
addition to requiring a non-linear minimization, a disadvan- 
tage of this method is that it requires evaluation over a large 
(3 dimensional) parameter space. 

The present invention utilizes a novel system that uses a 
chemometric analysis for measuring the fluorescence life- 35 
time and the fluorescence spectrum of a target mixture. 
Instead of irradiating the target sample with a single short 
pulse of light (photon counting) as other prior art systems, 
the present invention continuously irradiates the target mix- ^ 
ture with a light source whose amplitude (intensity) and 
phase (frequency) are modulated (varied) with time. In 
evaluating experimental phase resolved fluorescence spec- 
troscopy (PRFS) data, the variables sought are the 
concentrations, c-, the lifetimes, x i7 and the emission spectra, 45 
a,-(L), of the individual fluorescent contributors. The vari- 
ables that we may control include the emission wavelength, 

X, the modulation frequency, co, and the detector reference 
phase angle We first consider an amplitude modulated 
excitation: 

E(t)=Ej[l+M e sin ((Of)], *>0 (5) 

where co is the modulation angular frequency, M e is the 
modulation amplitude and E 0 is the average amplitude of the 55 
excitation. 

At steady state, the modulated part of the fluorescent 
intensity may be given by 

I K (t)=E a M < A sin((Of-<J>) (6) 60 

where A is the amplitude response at the fluorescent wave- 
length and T> is the phase lag at this wavelength and 
frequency. 

We originally detected this signal (the emitted signal) by 65 
correlating it with a reference signal (the driving/reference 
signal) by mixing the two signals with the mixing signal and 


integrating them over some multiple of a period. The cor- 
related signal could therefore be expressed as: 

a r% (7) 

S x (u), <b R ) = — / A (r)sin(wr - <b R )dt 

In Jo 

= E 0 M e Aco$$> R -<t>) 

1 v- 

= - E 0 M e 2_^ s i co*p i cos(fp i -& R ) 

1 i= l 

where r is the number of components present (or 
fluorophores) and is the phase of the reference signal 

relative to the backscattered excitation signal. We observe 
that S^cOjOj,) may be written as: 

S K (oy, <& r )=V 2 E J\f e [ G (co, A)co s <1^+5 (co, A) sin <!>*] (8) 

Therefore, as an alternative to Eqn. (8), we may determine 
G(co,X) and S(co,X) to within a common normalization factor, 
E 0 M e /2, for any fixed co and X by making measurements at 
two different, but known, values of In the presence of 
noise, more reference phase values would be used in the 
measurements and fit in a least squares manner to G and S. 

We then find that for each value of co and X (emission 
wavelength) we have (modulo the normalization factor): 


G{a), A) = ^ £ l (A)cos 2 9 ? j (ct>) 

(9) 

S(o), A) = ^ £ j (A)cos^ j (ct>)sin^ j (w) 

(10) 


with tan c^—cot:,-, where i runs over the number of fluoro- 
phores present and assumed to be acting in an independent 
(uncorrelated) fashion. 

The preferred embodiment of the present invention allows 
us to determine G and S by mixing the emitted signal from 
the mixture down to an intermediate frequency of 10 kHz, 
sample and digitize it at 50 kHz, Fourier transforming the 
result, and picking off the signal at 10 kHz. When this is 
done the mixed signal has the form: 

sin(w fF t)-S cos (co /F f)] (11) 

From which we see that G and S are proportional to the 
Fourier transform of Equation (11). 

In the case of noiseless data, where r fluorophores are 
present, it is necessary to make measurements at r values of 
L and r values of co to obtain: 

G=cc'A and S=cs-A, (12) 

where 

G =[ G jk\ 

s=[s jk \ 

CC=[COS 2 ^((Oy^fCCyJ 

cs=[cos ^-((Dyjsm <|),((0 / )]=[csy,] 

(13) 

where the quantities defined in boldface are matrices. The 
matrix indices i, j, and k run over the components, the 
modulation frequencies, and the emission wavelengths, 
respectively. 
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If no noise is present, we may take the indices i, j, and k 
above to run from 1 to r, the number of components. As the 
fluorophores are assumed to be independent and as the 
frequencies {oo y -} may be chosen arbitrarily, all the matrices 
defined above have true inverses. Thus, from Eqns. (12) and 
(13) we find: 


S=csA=wcc-TA 


w~ 1 SA~ 1 =cc m T=GA~ lm T 


S=cs-cc lm G 


S'G 1- CC=C5 


Now, we note that: 


cs = [cos^;(Wj)sin^;(toj) 


: [cos 2 (fi (w j)tan<pi (o j) 


By substitution from tan we have: 


10 G~ 1 -w~ 1 -SA~ 1 =A~ 1 -T (26) 

Solving this eigenequation determines t,- as the eigenval- 
ues and A -1 as the eigenvectors. However, as the eigenvec- 
tors are normalized, A may be better determined as above if 
!5 information about the component concentrations is required. 

Once the spectrum and lifetime of each component is 
known, we may use the information to identify the fluoro- 
phore. The concentration of each fluorophore may then be 
determined by comparing its spectrum against its standard 
20 spectrum. We find: 

c t (A) / ct i (A)t; (27) 

cj £j(X)/aj(X)Tj 


CS = [(cos <Pi{U)j)) 2 (x)jTi\ 




where w=diag[co ] and T=diag[xJ. Hence, we find 


S-G 1 -cc=wcc-T 


w ■ S-G 'CC=CC'T 


This equation is in the form of a standard eigenvalue/ 
eigenfunction equation, where the x t are the eigenvalues and 
the columns of cc are the eigenfunctions. We may also write 
this in the form: 


cc 1 -w lm S-G x -cc=T 


From which we have: 


det(w 1 ■ S ■ G 1 ) = r; 


30 where cq(X) is the standard spectral value of the ith com- 
ponent at wavelength \ and t,- is its lifetime. 

The number of components in a mixture may not be 
known a priori, and the method requires the number of 
frequencies and wavelengths employed to at least equal the 
35 number of components. Furthermore, some degree of noise 
will be present in the experimental data. Therefore, it is 
useful to determine a means of solving for lifetimes and 
spectra when the number of frequencies and the number of 
emission wavelengths are unequal and when each is greater 
40 than the number of independent components present in the 
mixture. 

In the following, we shall therefore assume that r com- 
ponents are present in a mixture, for which we have taken 
45 data at N frequencies and M emission wavelengths. We 
assume that N,M^r and that N^M, in general. 

From Eqns. (12) and (16), we therefore have: 


Tr(w~ l ■S-G~ l )=Y J T i 


(NxM)(Nxr)*(rxM) 

and 


Since (jj/cDy) is known once x t is known, we may determine 
cc or cs and solve the equation: 


G=ccA or S=csA 


for A with a standard LU decomposition, which gives us 


A=T~ lm ccpi m w~ 1 S (22) 

where ccpi is the pseudo inverse of cc. 

Another, equivalent, way of writing the eigenfunction 
equations proceeds from the observation that: 


( NxM) • (NxN) • (Nxr) (rxr) (rxM) 
where the number of rows and columns of each matrix is 
indicated below it. In the noise free case, all the matrices are 
of rank r. In the presence of noise, G and S will be of rank 
N or M, whichever is smaller. 

In the absence of noise, all matrices are of rank r. We 
manipulate Eqns. (28) and (29) slightly, to write: 


GA~ x =cc 


(23) 


w 1 -S=cc-TA 


(30) 
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where I is the rxr identity matrix. 

We may now make a singular value decomposition on the 
first of Eqns. (30): 

g=u-c 1 -v t (31) 

(NxM)- (Nxr) -(rxr) -(rxM) 

where C 2 is the diagonal matrix of singular values of G and 
U and V are each column orthonormal, i.e.: 

U T -U=V T -V=I. (32) 

From Eqns. (31) and (32), we have: 

C t =U T -G-V (33) 


deAP^. 

Therefore has an inverse. A similar argument shows 
that P T/ is also invertible. 

5 From Eqns. (30), (33) and (34), we have: 


C 2 =Pu'T'P v and C 1 =P L/ /P V . 


Therefore, we may write: 


C^C^Py^TPy 


We define: 

C 2 =u r -W~\s-V (34) 

where C 2 is not diagonal, in general, but is rxr. 

It may easily be shown, by a process that is similar to 
Gram-Schmidt orthogonalization, that while, in general, 
neither matrix is invertible, cc has a left inverse and A has 
a right inverse: 

cclrcc=I and A: Ari=I (35) 

From Eqns. (30) and (31), we may therefore write: 

cc=U-P IT and A=P v V r (36) 


Eqn. (42) is the eigenvalue equation for the matrix, 
Cr 1 ^, with eigenvalues equal to the diagonal elements of 
T, and eigenvectors proportional to the column vectors of 
PpT 1 . The spectrum, A, may be determined from Eqns. (13) 
and (30) since cc is determined once the lifetimes are 
known. 

A second method for solving Eqns. (30) is to use Eqn. (31) 
to write the pseudo inverse of G: 

Gpi=V-C x ~ x -U T (43) 

so that we have: 


P J j=C 1 V r Ari and P v =cclrU-C 1 . 

where P^ and P v are each rxr. 

We may interpret the first of Eqns. (36) in terms of its 
column vectors as the representation of the r vectors of cc in 
the orthonormal basis of the r vectors of U: 




where Cj is the jth column of cc, U,- is the ith column of U 
and P £ - is the corresponding element of P^. From this 
expression and the linear independence of the c and U £ , we 
see that P^ is invertible. Linear independence of the c ■■ 
requires that: 


G Gpi=U U r and Gpi G=V V T , 

where U-U 7 ’ is and NxN matrix and V-' V r is an MxM matrix, 
and each of these matrices is of rank r. 

From Eqns. (30) and (35), we may write: 


Substituting this into the first of Eqns. (30), we have: 
Gpi ■ w~ x -SAri= V ■ V r Ari ■ T 

From Eqns. (30) and (36) we have: 


W- X -S=UPrfTPyT 


45 From Eqn. (32), we therefore find: 


w - 1 -S=w- 1 S-V-V I . 


ajCLj = 0 only if aj = 0 for all j = 1, . . . , r. 
i,j= 1 


Substituting this into Eqn. (45), we therefore have: 

Gpi- w~ 1 -5 ■ V- V r A ri= V- V r AriT 


From Eqn. (36), we see that Eqn. (37) implies: 

Zw,=o. 

ij= 1 

From the linear independence of the U,, we must have 

A OT 

y PijOSj = 0, for all i = 1, ... ,r. 

7=1 

Thus, if the condition of Eqn (37) is to be satisfied, the 
system of Eqns. (39) must have only the trivial solution. This 
is equivalent to the requirement that 


Eqn. (47) is the eigenvalue equation for the matrix, 
Gpi-w^-S, with eigenvalues equal to the diagonal elements 
of T and eigenvectors proportional to the columns of 
V*V r Ari. The spectra, A, may be found as described above. 

When determining the lifetimes, we can renormalize Eqn. 
(47) by using 

pseudoinv(g) w -1 S evec=evec T (48) 

where evec=V-V :r -Ari for T-diag(x 1? . . . , x r ) 

For the case when noise is present, let us assume that r 
components are actually present, and let us take N<M. The 
experimentally determined matrices, Gx and Sx, are there- 
fore of rank N. We wish to determine the NxM matrices, G 
and S, of rank r that minimize some suitably defined error. 
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We shall define the Euclidean norm of a matrix, M, to be 
given by: 

||M || 2 = Tr(M T • M) = Tr(M ■M t ) = Y j Yj (49) 


The error we define will depend on our knowledge of the 
statistics of the experimental data. We shall assume that the 
errors at each data point are uncorrelated with those at other 
data points and that the standard deviations from the sample 
means are the same for all data points. This is tantamount to 
the assumption that the statistics are independent of the 
modulation frequency and the emission wavelength for the 
ranges of the variables explored. Previous experience with 
experimental data tends to bear this assumption out. 

Given these assumptions, it is reasonable to define the 
error, % 2 , in terms of the equally weighted distance of the 
experimental points from their analytical counterparts, i.e.: 

X 2 =||Gx-G|| 2 +||iv- 1 -5x->v- 1 -|| 2 (50) 

where G satisfies Eqn. (31) and where, from Eqns. (40) and 
(46), we see that S satisfies: 

w- 1 -S=U-C 2 -V r (51) 

Therefore, from Eqns. (31), (50) and (51), we must find U, 
V, C 1 and C 2 to minimize: 

y^=\\Gx- U-C 1 ■ V T \ | 2 +| |w _1 m Sx- U- C 2 - V 7 ]! 2 (52) 

Minimizing Eqn. (52) will yield a least squares fit for the 
data. 

For Eqn. (52) to be a minimum, its partial derivatives with 
respect to each free variable must vanish. Applying this 
requirement with respect to the variables of C l9 we find: 

C ± =U t -GxV (53) 

or, from Eqn. (31), 

G=UU t GxVV t =Q u GxQ v (54) 

where 

Q^u- U T and Qy=VV T (55) 

Proceeding similarly for the variables of C 2 , we find: 

C^U^w^-Sx-V (56) 

or 

w~ lm S=Q u'W -1 ' Sx'Q v (57) 

or 

=U-U T -w~ 1 -S-Vx-V T (58) 

Eqns. (54) and (57) determine C 1 and C 2 once U and V are 
determined. We now proceed to determine these variables. 

By substituting for G and w -1 S from Eqns. (54) and (57) 
into Eqn. (50), we find that minimizing that expression 
implies that we must maximize: 

yHQuGx-QvlMQ^w^-Sx-Q^ 2 (59) 
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^Trigg^+miw-'SXw-'sy] (60) 

Observing that: 

^ Qu'Qu=Qu> (61) 

and 

GvHSvT Qv'Qy^Qv* ( 62 ) 

io we may rewrite Eqn. (59) in the alternative forms: 

r ]=Tr(M v -U-U T ) (63) 

or 

15 ^TrQAfjV-V 7 } (64) 

where 

M v =Gx-QyGx T +{w- 1 -SxyQy{w~ 1 -Sx) T (65) 

20 

and 

M u =Gx r -Q u -Gx+{w- 1 -Sx) r -Q u iw- 1 -Sx) (66) 

We note that M v is an NxN matrix, M^is an MxM matrix, 
25 and both matrices are symmetric. 

If we assume that V is known, we may determine U by 
requiring that it be chosen to maximize Eqn. (63) subject to 
the constraint imposed by Eqn. (32). We may solve for U by 
30 means of Lagrange multipliers: 

Define: 

Tl'-n -Tr[Ly(U T U-W (67) 

35 where L v is a symmetric rxr matrix formed from the 
multipliers. 

We may now extremize r|' with U unconstrained to find: 

MyU=ULy (68) 

40 

Since it is symmetric, L v may be diagonalized by a 
similarity transform with an appropriate orthogonal matrix. 
Therefore, we may write: 

45 MyU=UODO~ x (69) 

where D is an rxr diagonal matrix and O is an rxr orthogonal 
matrix. 

We may then write: 

50 

MyUO=UOD (70) 

Substituting Eqn. (70) into Eqn. (63) we find: 

ri ^TriMy-U-O-O-^U 1 ) 

55 

r )=Tr[(UO)D(UO) T ] 

=Tr[D-(U-0)-(U-0 )] 

or, finally, 

60 

r\=Tr(D). (70) 

Thus, to maximize r| we choose U to maximize D as 
determined by Eqn. (70). We may easily do this by observing 
65 that Eqn. (70) has the form of an eigenvalue equation. The 
matrix D yields r of the N possible eignevalues of the matrix 
M v . We therefore solve the eigenvalue equation for M v : 


or 
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MyU=U-E (71) 

taking the largest r eigenvalues of E for D and, without loss 
of generality, choosing U to equal the corresponding eigen- 
vectors of U. 5 

A similar procedure may be used to determine V from 
Eqn. (64) if U is known using the relationship Mu*V=V*Lu 
to find L jj. 

FIG. 8 is a block diagram of the preferred embodiment of 
the present invention that uses the above foundation as a 10 
basis for a chemometric analysis of the emitted signal to 
determine the fluorescence spectrum and fluorescence life- 
times. The present invention first acquires the data from the 
data signal and the reference signal to determine the fre- 
quency and wavelength vectors 132. Next, the present 
invention assumes an arbitrary but fixed number of 
fluorophores, r, 134. We then determine the initial G and S 
matrices 136. 

As the frequencies in w are large, we renormalize 20 
||w _1 *Sx|| to equal ||Gx|| so as to avoid skewing the results in 
favor of the Gx data 138. This amounts to a rescaling of the 
units for the frequencies and the lifetimes. This rescaling of 
units is compensated once the lifetimes are found, so they 

9,5 

are expressed in seconds. 

We begin the iterative approach by determining the initial 
approximations for U, V, g and w -1 s approximations 140, 
142, and 144. The initial iteration begins by choosing the r 
columns of U 1 to be the eigenvectors corresponding to the 30 


r largest eigenvalues of the matrix 


Gx-Gx r +{w~ 1 -Sxy {w- x -Sx) T 

(73) 

or 



(74) 


where L=diag(X 1 , . . . , X r ) forX 1 > . . . >\ r 
Next, U 1 is used to form the matrix M U1 in accordance with 
Eqn. (66). We can determine W 1 from M m by solving the 40 
eigenfunction equation M t/1 *V=V*L t/ by choosing its r col- 
umns to be the eigenvectors corresponding to the r largest 
eigenvalues of M^. We then form M V1 in accordance with 
Eqn. (65). 

We next initialize the decision variables 146 to produce 45 
the solution in a reasonable time, with a reasonable accuracy, 
and a reasonable number of iterations. We begin the iterative 
solution by comparing the decision variables 148 for con- 
vergence to see if further determinations for U, V, g and w -1 s 
are necessary as above. In the iterative solution, we proceed 50 
to maximize r\, given by Eqn. (59), and recognizing that: 

0^ri=l|Gx|| 2 +||iv- 1 -5x|| 2 =JV (72) 

Proceeding in this way, we form a sequence of r\ values 55 
{r|,-}. We see easily that the sequence is mono tonic non- 
decreasing. By Eqn. (72) and the Bolzano-Weierstrass 
theorem, the sequence has at least one limit point in the 
interval [0,N]. The mono tonic nature of the sequence guar- 
antees that there can be but one limit point. Therefore the 60 
sequence must converge. 

The values of U and V corresponding to the limit of the 
sequence are then used to form G and w -1 S in accordance 
with Eqns. (54) and (57). Eqn. (47) may subsequently be 65 
solved for the lifetimes and spectra as described previously 

160, 162, and 164. 
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In the preceding, we assumed an arbitrary but fixed value, 
r, for the number of fluorescing components present. Fits to 
the data may be found for different values of r. The question 
then arises as to the value of r that best fits the experimental 
data. Two qualitative criteria may be stated: (1) the spectra 
should look reasonable; and (2) the lifetimes should be 
positive real numbers. 

Somewhat more quantitatively, we recognize that % 2 , 
suitably renormalized for frequency as described above, 
represents our error for a fixed value of r. However, the value 
of x 2 depends on the number of free parameters available for 
the fit, which, in turn, depends on r. Instead, we attempt to 
remove this dependence by dividing by the difference 
between the number of experimental data points and the 
number of free parameters present: 


where y? R is the reduced error term, N^ is the number of 
experimental data points and N f is the number of free 
parameters. 

Since we collect an amplitude and a phase value for each 
emission wavelength and each modulation frequency, we 
have: 


Nx=2N^N x (76) 

where Nw is the number of frequencies and N x is the number 
of wavelengths used. 

From Eqns. (54), (57) and (68), we see that the free 
parameters being fit to the experimental data are the ele- 
ments of U, V, L v and L u subject to the constraints of Eqns. 
(32). For r components, we therefore find that: 

1 1 (77) 

N f = r- Nw + r-N x +2- - -r(r+ l)-2- --r(r+ 1) 

= r ■ N w + r-Nx 


We evaluate % 2 R for each r, and choose the value of r 
corresponding to the smallest value of % 2 R to represent the 
number of components which best fits the data 158. 

If two or more lifetimes are the same, their components 
will not be separable by this method, nor will they be 
separable by the method of PREEMs. Their contributions 
will remain combined as that of a single component with the 
common lifetime and a spectrum given by the combined 
spectrum appropriate to their concentrations. 

FIG. 9 illustrates an embodiment of the present invention 
for the method described in FIG. 8 for a chemometric 
analysis of the emitted signal to determine the fluorescence 
spectra and fluorescence lifetimes. The user interface 170 
allows the user of the present invention to interact with the 
system of the present invention. The user interface allows 
the user to choose the number of components or fluoro- 
phores to solve for. The data acquisition component 172 uses 
the processor 50 coupled to a data gathering component 
through a data input/output path 68 to acquire the data from 
the emitted signal to determine the frequency and wave- 
length vectors. A determiner 174 determines the G and S 
matrices from the frequency and wavelength vectors. A 
renormalizer 176 then renormalizes w _1 S^. to G^. A deter- 
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miner 178 next determines the initial U approximation, and 
a determiner 180 determines the initial V approximation, 
and a determiner 182 determines the initial g and w -1 s 
approximations. An initializer/updater 183 then initializes 
the decision variables. The iterative solution of the present 5 
invention uses a comparator 184 that compares the decision 
variables for convergence to see if further determinations for 
U, V, g and w -1 s are necessary. If the solution does not 
converge, the present invention uses the determiner 178 to 10 
determine the next U approximation, the determiner 180 to 
determine the next V approximation, and the determiner 182 
to determine the next g and w -1 s approximations. The 
present invention then uses the initializer/updater 183 to 
update the decision variables. 15 

When the solution converges, a determiner 186 deter- 
mines the reduced best fit error for the analysis of the 
individual fluorescence lifetimes and fluorescence spectra. 
The present invention then uses an extractor 188 that 20 
extracts the individual fluorescence lifetimes from the target 
mixture, and an extractor 190 that extracts the individual 
fluorescence spectra. 

FIG. 10 illustrates the present invention’s ability to iden- 
tify and discriminate between individual overlapping spec- 25 
tral components in a target mixture. The present invention 
allows for the extraction of both the individual fluorescence 
spectra and lifetimes from a target mixture of fluorophores. 
FIG. 10 illustrates the present inventions ability to differ- 3Q 
entiate spectra using a three-dye mixture of OXAZINE™ 
720, 725 and 750. OXAZINE is a trademark of Excition, 
Inc. We mixed the dyes at equal concentrations of 3.3 //M. 

We recorded the emission spectra for an excitation signal 
wavelength of 640 nm and the laser modulation frequency 35 
was swept (modulated) from 10 MHz to 140 MHz at 5 MHz 
increments. The present invention extracted fluorescence 
lifetimes of 3.705 nsec for Oxazine 720 (196 on FIG. 10 ), 
1.979 nsec for Oxazine 750 ( (198 on FIG. 10 ), and 0.5588 
nsec for Oxazine 725 (200 on FIG. 10 ). The results from the 
present invention compared very well with the listed indi- 
vidual dye fluorescence lifetimes of 3.739 nsec for Oxazine 
720, 2.014 nsec for Oxazine 750 and 0.9935 nsec for 
Oxazine 725. The individual spectra extracted for each dye 45 
from the mixture revealed spectral characteristics that 
matched with spectra obtained from the individual dyes. 

The present invention overcomes the limitations of the 
prior art systems by utilizing a novel technique to measure 
the fluorescence lifetime. Instead of irradiating the target 
sample with a single short pulse of light, the present inven- 
tion continuously irradiates the target sample with a light 
source whose amplitude and phase are modulated with time. 
This technique allows the present invention to use a chemo- 55 
metric analysis to automatically extract the lifetimes from 
the ‘phase delay’ and ‘intensity vs. time’ characteristics of 
the emitted light. 

The present invention is a system for chemometric analy- 
sis for the extraction of the individual fluorescence spectrum 60 
and fluorescence lifetime from a target mixture. The present 
invention comprises a processor with an apparatus for 
generating an excitation signal to transmit at a target mixture 
and an apparatus for detecting the emitted signal from the 65 
target mixture. The present invention extracts the individual 
fluorescence spectrum and fluorescence lifetime measure- 
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ments from the frequency and wavelength data acquired 
from the emitted signal. The present invention uses an 
iterative solution that first requires the initialization of 
several decision variables and the initial approximation 
determinations of intermediate matrices. The iterative solu- 
tion compares the decision variables for convergence to see 
if further approximation determinations are necessary. If the 
solution converges, the present invention then determines 
the reduced best fit error for the analysis of the individual 
fluorescence lifetime and the fluorescence spectrum before 
extracting the individual fluorescence lifetime and fluores- 
cence spectrum from the emitted signal of the target mixture. 

The present invention additionally comprises a method 
and apparatus for generating and mixing signals for 
frequency-domain lifetime and spectral fluorometry. The 
present invention comprises a plurality of signal generators 
that generate a plurality of signals where the signal genera- 
tors modulate the amplitude and/or the frequency of the 
signals. The present invention uses one of these signals to 
drive an excitation signal that the present invention then 
directs and transmits at a target mixture, which absorbs the 
energy from the excitation signal. The property of fluores- 
cence causes the target mixture to emit an emitted signal that 
the present invention detects with a signal detector. The 
present invention uses a plurality of mixers to produce a 
processor reference signal and a data signal. The present 
invention then uses a processor to compare the processor 
reference signal with the data signal by analyzing the 
differences in the phase and the differences in the amplitude 
between the two signals. The processor then extracts the 
fluorescence lifetime and fluorescence spectrum of the emit- 
ted signal from the phase and amplitude information using 
a chemometric analysis. 

Other embodiments of the invention will be apparent to 
those skilled in the art after considering this specification or 
practicing the disclosed invention. The specification and 
examples above are exemplary only, with the true scope of 
the invention being indicated by the following claims. 

We claim: 

1. A method of chemometric analysis for the extraction of 
an individual fluorescence spectrum and an individual fluo- 
rescence lifetime from a target mixture comprising the 
following steps: 

acquiring the frequency and wavelength data from the 
emitted signal from the target mixture; 

determining a G matrix and a S matrix from a frequency 
vector and a wavelength vector; 

renormalizing w -1 S^ to G^; 

determining an initial U, V, g and w -1 s approximations; 

initializing one or more decision variables; 

comparing the decision variables for convergence to see 
if further determinations for U, V, g and w -1 s are 
necessary; 

extracting the individual fluorescence lifetime from the 
target mixture; and 

extracting the individual fluorescence spectrum from the 
target mixture. 

2. The method of claim 1 wherein the step of determining 
the initial U, V, g and w -1 s approximations further com- 
prises the following steps: 
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determining the initial U approximation; 
determining the initial V approximation; and 
determining the initial g and w -1 s approximations. 

3. The method of claim 1 wherein the step of comparing 
the decision variables for convergence further comprises the 5 
following steps: 

determining a next U approximation; 

determining a next V approximation; 

determining a next g and w -1 s approximations; and 10 

updating the decision variables. 

4. The method of claim 1 further comprising the step of 
determining a best fit error for an analysis of the individual 
fluorescence lifetime and the fluorescence spectrum. 

5. A program storage device readable by a computer, 
tangibly embodying a program of instructions executable by 
the computer to perform method steps for a method of 
chemometric analysis for the extraction of an individual 
fluorescence spectrum and an individual fluorescence life- 20 
time from a target mixture, comprising the following method 
steps: 

acquiring the frequency and wavelength data from the 
emitted signal from the target mixture; 
determining a G matrix and a S matrix from a frequency 25 
vector and a wavelength vector; 
renormalizing w -1 S x to G x ; 

determining an initial U, V, g and w -1 s approximations; 
initializing one or more decision variables; 30 

comparing the decision variables for convergence to see 
if further determinations for U, V, g and w -1 s are 
necessary; 

extracting the individual fluorescence lifetime from the 
mixture; and 35 

extracting the individual fluorescence spectrum from the 
target mixture. 

6. The program storage device of claim 5 wherein said 
step of determining the initial U, V, g and w -1 s approxima- ^ 
tions further comprise the following steps: 

determining the initial U approximation; 
determining the initial V approximation; and 
determining the initial g and w -1 s approximations. 

7. The program storage device of claim 5 wherein said 45 
step of comparing the decision variables for convergence 
further comprises the following steps: 

determining a next U approximation; 
determining a next V approximation; 50 

determining a next g and w -1 s approximations; and 
updating the decision variables. 

8. The program storage device of claim 5 further com- 
prising the step of determining a best fit error for an analysis 

of the individual fluorescence lifetime and the fluorescence 55 
spectrum. 

9. A system for chemometric analysis for the extraction of 
an individual fluorescence spectrum and an individual fluo- 
rescence lifetime from a target mixture, comprising: 

a processor that acquires the frequency and wavelength 
data from the emitted signal from the mixture; 
means for determining a G matrix and a S matrix from a 
frequency vector and a wavelength vector; 
means for renormalizing w _1 S JC to G^; 65 

means for determining an initial U, V, g and w -1 s approxi- 
mations; 
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means for initializing one or more decision variables; 

means for comparing the decision variables for conver- 
gence to see if further determinations for U, V, g and 
w -1 s are necessary; 

means for extracting the individual fluorescence lifetime 
from the target mixture; and 

means for extracting the individual fluorescence spectrum 
from the target mixture. 

10. The system of claim 9 wherein said means for 
determining the initial U, V, g and w -1 s approximations 
further comprise: 

means for determining the initial U approximation; 

means for determining the initial V approximation; and 

means for determining the initial g and w -1 s approxima- 
tions. 

11. The system of claim 9 wherein said means for 
comparing the decision variables for convergence further 
comprises: 

means for determining a next U approximation; 

means for determining a next V approximation; 

means for determining a next g and w -1 s approximations; 
and 

means for updating the decision variables. 

12. The system of claim 9 further comprising means for 
determining a best fit error for an analysis of the individual 
fluorescence lifetime and the fluorescence spectrum. 

13. An apparatus for chemometric analysis for the extrac- 
tion of an individual fluorescence spectrum and an indi- 
vidual fluorescence lifetime from a target mixture, compris- 
ing: 

a processor that acquires the frequency and wavelength 
data from the emitted signal from the mixture; 

a determiner that determines a G matrix and a S matrix 
from a frequency vector and wavelength vector; 

a renormalizer that renormalizes w -1 S^ to G^; 

a determiner that determines an initial U, V, g and w -1 s 
approximations; 

an initializer that initializes one or more decision vari- 
ables; 

a comparator that compares the decision variables for 
convergence to see if further determinations for U, V, g 
and w -1 s are necessary; 

an extractor that extracts the individual fluorescence life- 
time from the target mixture; and 

an extractor that extracts the individual fluorescence spec- 
trum from the target mixture. 

14. The method of claim 13 wherein said determiner that 
determines the initial U, V, g and w -1 s approximations 
further comprises: 

a determiner that determines the initial U approximation; 

a determiner that determines the initial V approximation; 
and 

a determiner that determines initial g and w -1 s approxi- 
mations. 
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15. The method of claim 13 wherein said comparator that 
compares the decision variables for convergence to see if 
further determinations are necessary further comprises: 
a determiner that determines a next U approximation; 
a determiner that determines a next V approximation; 
a determiner that determines a next g and w -1 s approxi- 
mations; and 
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an updater that updates the decision variables. 

16. The method of claim 13 further comprising a deter- 
miner that determines a best fit error for an analysis of the 
individual fluorescence lifetime and the fluorescence spec- 
trum. 



